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We report results on the branching fraction (BJ 7 ) measurement of the color-suppressed 
decays B° -> DV, D* tt , D°r], D*°r], D°lu, D*°u, D°r}', and D*°rf. We measure the 
branching fractions BJ r (D°n°) = (2.78 ± 0.08 ± 0.20) x 10" 4 , 
BJ~(D*°iT ) = (1.78 ± 0.13 ± 0.23) x 10" 4 , BJ z (D Q rj) = (2.41 ± 0.09 ± 0.17) x 10~ 4 , 
BFiD* !]) = (2.32 ± 0.13 ± 0.22) x 10" 4 , BF(D°lo) = (2.77 ± 0.13 ± 0.22) x 10" 4 , 
BF{D*°u) = (4.44 ± 0.23 ± 0.61) x 10" 4 , BT{D°i) = (1.38 ± 0.12 ± 0.22) x 10~ 4 and 
BJ r (D*°r]') = (1.29 ± 0.23 ± 0.23) x 10" 4 , where the first uncertainty is statistical and the 
second is systematic. The result is based on a sample of (454 ± 5) x 10 6 BB pairs collected 
at the T(4S) resonance from 1999 to 2007, with the BABAR detector at the PEP-II storage 
rings at the Stanford Linear Accelerator Center. The measurements are compared to 
theoretical predictions by factorization, SCET and pQCD. The presence of final state 
interactions is confirmed and the measurements seem to be more in favor of SCET 

compared to pQCD. 
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1 INTRODUCTION 



Weak decays of hadrons provide a straight access to the parameters of the CKM matrix and thus 
to the study of the CP violation. Gluon scattering in the final state (Final State Interactions, 
or FSI) can modify the decay dynamics and so must be well understood. The two-body hadronic 
decays with a charmed final state, B — > Dh, are of great help in studying strong-interaction physics 
related with the confinement of quarks and gluons into hadrons. 

The decays B — > Dh, where h is a light meson, can proceed through the emission of a 
boson following three possible diagrams: external, internal (see Fig. [I]) or by a boson exchange 
whose contribution is negligible pp. 




d * d 




d « d 



Figure 1: External (a) and internal (b) tree diagrams for B° — ► Dh decays. 

In the case of the decays B° -» D^°h°, the major contribution comes from the internal 
diagram [2]. Since mesons are color single objects, in internal diagrams 

B o D (*)0 h the quar k s 

from the decay are constrained to have the anti-color of the spectator quark, which induces a 
suppression of internal diagrams in comparison with external ones. Internal diagrams are so called 
color-suppressed and external ones are called color- favored. 

In the factorization model [21 S HI [5] , the non-factorizable interactions in the final state by soft 
gluons are neglected. The matrix element in the effective weak Hamiltonian of the decay B — > Dh 
is then factorized into a product of asymptotic states. Factorization appears to be successful in the 
description of the color-favored decays [6] . 

The decays B° -» D^° 7T° were first observed by CLEO [7] and Belle [8] with respectively 9.67 
and 23.1xl0 6 BB pairs. The Belle collaboration has also observed the decays D°rj and D°uj and 
put upper limits on the BJ- of D*°rj and D*°lo [8]. 

The branching fraction (BJ-) of the color-suppressed decays B° — > D^°tt°, D^°r], D^*'°u, and 
D°r)' were measured recently by BABAR [9\ with 88 x 10 6 BB pairs and an upper limit was set on the 
BT of D*°r]'. The Belle collaboration measured with 152 x 10 6 BB pairs the BT of B° -» D^°h°, 
h° = 7T°, 77, u, and rj' pU [IT] and studied the decays B° -> D^°p° with 388 xlO 6 BB pairs [12] . 
Many of these measurements showed a significant disagreement with predictions by factoriza- 
tion [13J, but stronger experimental constraints are needed to distinguish between the different 
models of the color-suppressed dynamics like pQCD (perturbative QCD) [HI [15] or SCET (Soft 
Collinear Effective Theory) |16 (. I17 |, [T5] . This paper reports the branching fraction measurement of 
eight color-suppressed decays B° -» L>W 7r°, D^°r], Z?(*)°(j and D^°rf with 454xl0 6 BB pairs. 



8 



2 THE BABAR DETECTOR AND DATASET 



The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e + e~ storage ring. The BABAR detector is described in detail in Ref. [19] . Charged particle tracks 
are reconstructed using a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) 
immersed in a 1.5 T magnetic field. Tracks are identified as pions or kaons (particle identification or 
PID) based on likelihoods constructed from energy loss measurements in the SVT and the DCH and 
from Cherenkov radiation angles measured in the detector of Cherenkov light (DIRC). Photons are 
reconstructed from showers measured in the electromagnetic calorimeter (EMC). Muon and neutral 
hadron identification are performed with the instrumented flux return (IFR). 

The data sample consists of an integrated luminosity of 413 fb _1 recorded at the T(4S) reso- 
nance with a center-of-mass (CM) energy of 10.58 GeV, corresponding to (454 ±5) x 10 6 BB pairs. 
A data sample of 41 fb™ 1 with a CM energy 40 MeV below the T(4S) resonance is used to study 
background contributions from continuum events e + e~ — > qq (q = u, d, s, c). 

Samples of simulated Monte Carlo (MC) events were used to determine signal and background 
characteristics, optimize selection criteria and evaluate efficiencies. Simulated events e + e~ — > 
T(AS) — > B + B~ , B°B°, e + e~ — > qq (q = u, d, s) and e + e~ — > cc are generated with EvtGen |20j, 
which interfaces to Pythia [2T] and Jetset [22]. Separate samples of exclusive B° — ► D^*'°h° 
decays were generated to evaluate the signal features and the efficiency of selections on signal. A 
sample of exclusive B~ — > was generated for the study of that background. All MC samples 

include simulation of the BABAR detector response generated through GEANT4 [23j. The integrated 
luminosity of the MC samples is about three times the data luminosity for BB, one times the data 
luminosity for e + e~ — > qq (q = u, d, s) and two times for e + e~ — > cc. The equivalent integrated 
luminosities of the exclusive simulations range from 50 to 2500 times the data luminosity. 

3 ANALYSIS METHOD 

3.1 Event reconstruction 

Charged particles tracks are reconstructed from measurements in the SVT and/or the DCH, and 
an identification is assigned by the PID algorithm. Extrapolated tracks must be in the vicinity of 
the e + e~ interaction point, i.e. within 1.5 cm in the plane transverse to the beam axis and 2.5 cm 
along the beam axis. The tracks used for the reconstruction of rj — ► 7r + 7r~7r must in addition have 
a transverse momentum px larger than 50 MeV/c. When PID criterion is required on a track, the 
track polar angle 9 must be in the DIRC fiducial region 25.78° < 9 < 146.10°. Photons are defined 
as single bumps in the EMC crystals not matched with any track, and with a shower lateral shape 
consistent with a photon. Because of high background in the forward region of the EMC caused 
by the beam asymmetry, the photons detected in the region 9 < 21.19° are rejected. 

Intermediate resonances of the decays B° — ► D^*'°h° are reconstructed by combining tracks 
and/or photons for the channels with the highest decay rate and detection efficiency. When mesons 
are combined to build a resonance, masses are fixed to their nominal value [24] . For the u and 
p° mesons, whose natural width is not negligible in comparison with the experimental resolution, 
their mass is not fixed to the nominal value. The selections applied to each meson ir°, rj, uj, t/, 
D°, and D*° are optimized by maximizing the figure of merit S/^J S + B where S is the number 
of signal and B is the number of background events. The numbers S and B are computed from 
simulations, the BJ-'s used to evaluate S are the world average value given by PDG [23]. Each 
resonance mass distribution is fitted with a set of Gaussian functions or a modified Novosibirsk 
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function [29], which is composed of a Gaussian-like peaking part with two tails at low and high 
values. Resonance candidates are then required to have a mass within ±2.5 a around the fitted 
mass central value, where a is the resolution of the mass distribution obtained by the fit. For the 
resonances D° - K- ~it + it and D*° — > D°j, the lower bound is extended to —3a because of the 
photon energy losses in front and between the EMC crystals, which makes the mass distribution 
asymmetric with a tail at low values. 

3.2 Selection of intermediate resonances 

3.2.1 7T° selection 

The 7T° mesons are reconstructed by combining two photons, each photon energy £'(7) must be 
larger than 85 MeV for 7r° coming from B° decay, and larger than 60 MeV for ir° coming from 77, 
uj or D°. Soft 7r°'s coming from D*° — ► D°tt° must satisfy £'(7) > 30 MeV. The reconstruction 
resolution for high momentum 7T° iS limited by the angle between the two daughter photons; for 
low momentum 7r° the resolution is limited by the neutral hadron background in EMC. The tt° 
reconstructed mass resolution is about 6 MeV/c 2 for tt° coming from 77, 7 MeV/c 2 for 7r° coming 
from uj or D°, and 8 MeV/c 2 for tt° coming from B° or D*°. The selection efficiency on signal 
ranges from 85 to 93 %. 

3.2.2 77 selection 

The 77 mesons are reconstructed in the 77 and tt + tt^tt° decay modes. These modes account for 
62 % of the total decay rate [24] , and may originate from B° — > D^*' G r] or 77' —* ir + ir~r] decays. 

The rj — > 77 candidates are reconstructed by combining two photons that satisfy E{^f) > 
200 MeV for B° daughters and E{^) > 180 MeV for 77' daughters. As high momentum tt°s may 
fake 77 — > 77 decays, a veto is applied against ir°: for each 77 — > 7172 candidate, the photons 71/2 
are associated with photons of the rest of event 7^. If E{^i) > 200 MeV and the invariant mass 
of the pair {7171} or {727^} is in the ir° mass window 115 < ?™(7i/27i) < 150 MeV/c 2 , then the 77 
candidate is rejected. The resolution of the 77 — > 77 mass distribution is dominated by the EMC 
resolution and is about 15 MeV/c 2 . 

For 77 candidates reconstructed in 7r + 7r~7r°, the ir° is required to satisfy the conditions described 
in Section \3. 2. 11 The mass resolution is about 3.5 MeV/c 2 , which is smaller than for 77 — > 77 thanks 
to the higher resolution of the tracking system. The selection efficiency on signal is about 77 % for 
77 — > 77 and 75 % for 77 — > -k + tt~it . 

3.2.3 u selection 

The uj mesons are reconstructed in the 7r + 7r _ 7r° decay mode. These modes account for 89 % of the 
total decay rate. The ir°s are required to satisfy the conditions described in Section 13.2.11 and the 
must fulfill the condition pt^^ 1 ) > 200 MeV/c. The natural width of the uj mass distribution 
r ~ 8.49 MeV |24j is comparable to the experimental resolution a ~ 7 MeV/c 2 , therefore the uj 
mass is not constrained to its nominal value. We define a total width a tot = V^ 2 + T 2 ~ 11 MeV/c 2 
and require the uj candidates to satisfy \m(uj) — m(u;)mean| < 2.5 otot- The selection efficiency on 
signal is about 82 %. 
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3.2.4 p° selection 

The p° mesons originate from rj — ► p 0r y and are reconstructed in the 7r + 7r _ decay mode. These 
modes account for 100 % of the total decay rate. Charged particle tracks must satisfy Pt{^) > 
100 MeV/c. We define the helicity angle 9 p o as the angle between the pion momentum in the p° 
rest frame and the p° momentum in the rj rest frame. Because the p° is a vector and pion is a 
pseudoscalar, the angular distribution is proportional to sin(# p o) 2 for pure signal and is flat for 
background. The p° candidates with |cos(# p o)| > 0.73 are rejected. Due to the p° large natural 
width r ~ 149.4 MeV [23], no mass constraint is applied to the p°. The mass of the p° candidates 
must be within 160 MeV/ c 2 around the nominal mass value. 

3.2.5 rj selection 

The 7/ mesons are reconstructed in the tt + tt^t](^ 77) and p 0r y decay modes. These modes account 
for 30 % of the total decay rate. For the reconstruction of rj — > n + tt~ ri(-^>- 77), the rj candidates 
must satisfy the selections described in Section [3.2.21 The rj mass resolution is about 3 MeV/c 2 . 

For rj candidates reconstructed in p°7, the /j 0, s are required to satisfy the conditions described 
in Section 13.2.41 and the photons must have an energy larger than 200 MeV. As photons coming 
from 7T° decays may fake signal, the veto against 7r° described in Section 13.2.21 is applied. The rj 
mass resolution is about 8 MeV/c 2 , which is worse than rj —* 7r + 7r~7/ because of the resolution on 
gamma reconstruction and the large p mass width. 

The selection efficiency on signal is about 69 % for rj — > Tr + ir~r}(— > 77) and 66 % for rj — > p°7. 

3.2.6 K° s selection 

The Kg mesons are reconstructed in the ir~ir + decay modes. These modes account for 69 % of 
the total decay rate. The x 2 probability of the vertex fit of charged pions must be larger than 
0.1 %. We define the flight significance as the ratio L/cjl where L is the K® flight length in the 
plane transverse to the beam axis and gl is the uncertainty on L determined from the vertex fit 
constraint. The combinatorial background is rejected by requiring a flight significance larger than 
5. The reconstructed K® mass resolution is about 2 MeV/c 2 . The selection efficiency on signal is 
about 86 %. 

3.2.7 D° selection 

The D° mesons are reconstructed in K~n + , K~Tr + ir°, K^tt + tt^tt + , and K®ir + Tr~ decay modes. 
These modes account for about 28 % of the total decay rate. All D° candidates must satisfy 
p*(D°) > 1.1 GeV/c. That requirement is loose enough that background populate the sidebands 
of the signal region. The tt^'s coming from the D° candidate must fulfill pr{^) > 400 MeV/c 
for D° -> K~7T + , pt(vt ± ) > 100 MeV/c for D° -> K-tt + tt-tt + and pt(vt ± ) > 120 MeV/c for 
D° — > K®7r~7r + . The x 2 probability of the vertex fit of charged pions must be larger than 0.1 % for 
D° — > K~ir + and larger than 0.5 % for the other modes. The kaon candidates must satisfy tight 
kaon criteria for K~tt + and K^ir + ir , and a loose kaon criteria for K~ir + ir~7r + . For KgTT + ir~ , the 
Kg candidates must satisfy the selection criteria described in Section 13.2.61 

The decay D° to K~ ~7r + 7T° proceeds mainly through the resonances K* (K*° — > K ir + or K* — > 
K~ir°) and p + (7r + 7r°). Combinatorial background is rejected by using the parametrization of the 
K~tt + tt° Dalitz distribution studied by the Fermilab E691 experiment [25J. The D° candidates 
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that are not in the resonance regions of the Dalitz distribution are rejected. The ir° must satisfy 
the selections described in Section 13,2.11 

The reconstructed D° mass resolution is about 4, 5, 6, and 11 MeV/c 2 for K^tt + tt^tt + , Kg-K + ir^ , 
K~tt + , and K~ir + ir° modes, respectively. To account for the asymmetry in the mass distribution of 
~ir + ir°, we require that the D candidates satisfy —3 a < (m(D°) — rn(D°) mean ) < 2.5 a. 

The selection efficiency on signal is about 71 % for K~n + , 60 % for K~-K + Tr~ir + , 71 % for 
K°tt + tt-, and 44 % for K--k + tt°. 

3.2.8 D*° selection 

The D*° mesons are reconstructed in the D°tt° and D°7 decay modes. The tt° and D° candidates are 
required to satisfy the selections described in Section [3. 2 . 1 1 and 13. 2 . 71 respectively. The photons from 
D*° — >■ D°7 must fulfill £'(7) > 200 MeV and pass the veto against ir° described in Section 13.2.21 

The resolution of the mass difference Am. = m(D*°) — m(D°) is about 2 MeV/c 2 for D°ir° 
and 7 MeV/c 2 for D°^f. The D*° — ► D°7r° candidates must satisfy |Am — Ammeanl < 2.5 <r. 
Because of the asymmetry of the Am distribution for D°7, the D*° candidates must satisfy —3 a < 
(Am - Ammean) < 2.5 a. 

The selection efficiency on signal is about 49 % for D°ir° and 39 % for Z?°7. 

3.3 Selection of B candidates 

The B candidates are reconstructed by combining a with an hP, with the D*° and h° masses 

constrained to their nominal value except when h° is the u. One needs to discriminate real B's 
from fake ones created from combinatorial or crossfeed background. 

3.3.1 B kinematic variables 

Two kinematic variables are used in BABAR to select B candidates: the energy-substituted mass twes 
and the energy difference AE. These two variables use the constraints from the precise knowledge 
of the beams' energies and from energy conservation in the two-body decay T(45) — * BB. The 
quantity mss is the invariant mass of the B candidate where the B energy is set to the beam energy 
in the CM frame: 



and AE is the energy difference between the reconstructed B energy and the beam energy in the 
CM frame: 



where y/s is the e + e~ CM energy and (E^,Pb) is the B quadrivector in the CM frame. For B signal 
events, the mES distribution peaks at the B mass with a resolution of about 3 MeV/ c 2 dominated by 
the beam energy spread, whereas AE peaks near zero with a resolution of 15 — 50 MeV depending 
on the number of photons in the final state. 




(1) 




(2) 
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3.3.2 Rejection of e + e — * qq background 

The continuum background e + e~ — > qq (q = u, d, s or mainly c) creates real high momentum 
mesons D^°, tt°, t/'\ uj that can fake the signal mesons originating from the two body decays 
B° — > D^*'°h°. That background is thus not rejected by the selections on intermediate resonances. 
Since the B mesons are produced almost at rest in the T(4S) frame, so the T(4f>) — * BB events 
shape is spherical. By comparison, the qq events have a back-to-back jet-like shape given that the 
quark masses (q 7^ b) are small compared to the CM energy. The qq background can hence be 
discriminated by the event shape described by the following variables: 



The thrust angle 9t defined as the angle between the thrust axis of the B candidate and the 
thrust axis of the rest of event. The thrust axis T is defined as the axis that maximizes the 
quantity T: 

T = E ij?J, r| ■ (3) 
The distribution of | cos(#t)| is flat for signal and peaks at 1 for continuum background. 
Legendre monomials Lq and L2 defined as: 

Lo = $>* ; L 2 = 5>*|cos(0*)| 2 , (4) 



with p* the momentum of the particle i that does not come from a B candidate, and 9* is 
the angle between p* and the thrust axis of the B candidate. 

• The polar angle 6* B between the B momentum in the T(4S) frame and the beam axis. The 
T(4S) being vector (J = 1) and the B mesons being pseudoscalar (J = 0), the angular 
distribution is proportional to sin 2 (#|j) for signal and roughly flat for background. 

These four variables are combined in a Fisher discriminant built with the TMVA [26] toolkit 
package. The Fisher J r s hape is trained with signal MC events and off-peak data events; in order to 
maximize the number of off-peak events all the B° — > D^°h° modes are combined. The training 
and test are performed with 2 x 20000 signal events and 2 x 20000 off-peak events. The obtained 
Fisher formula is: 

T BhaV e = 2.36-1.18- I cos(0 T )| + 
0.20 • L - 1.01 • L 2 - 0.80 • I cos(0|j)|. (5) 

The qq background is rejected by applying a selection on Fshape- The selection is optimized for 
each B signal mode by maximizing the statistical significance with signal MC against generic MC 
e + e~ — ► qq, q 7^ b. That requirement retains between 36 % and 98 % of B signal, while rejecting 
between 23% and 97 % of the qq background respectively. 

3.3.3 Rejection of T(45) -» BB background 

The u mesons in B° — > D°lu are polarized. We define the normal angle On [HI E2] as the angle 
between the normal to the to decay plane and the direction of B° in the to rest frame. That definition 
is the equivalent of the two-body helicity angle for the three-body decay. To describe the 3-body 
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decay distribution of uj — ► 7r + 7r~7r°, we define the Dalitz angle Op [9j as the angle between the ir° 
momentum in the u frame and the tt + momentum in the {7r + 7r - } frame. 

Given the angular momenta of the B° and D° mesons (J = 0), and uj meson (J = 1), the 
signal distribution is proportional to cos 2 (#jv) and sin 2 (#£)) while the combinatorial background 
distribution is roughly flat. These two angles are combined in a Fisher discriminant Thel built from 
signal MC events and generic qq and BB MC events: 

T hel = -1.41 - 1.01 • | cos(0£>)| + 3.03 • | cos(0jv)|. (6) 

We require B° — > D°u> candidates to satisfy Thel > —0.1, to obtain an efficiency on signal (back- 
ground) of 85 % (39 %). 

We also use the angular distribution of D*° — > D°7r° to reject combinatorial background. We 
define the helicity angle 9d* as the angle between the D° momentum in the D*° frame and the 
D*° momentum in the B° frame. The angular distribution is proportional to cos(#£)*) 2 for signal 
and roughly flat for combinatorial background. Selection on |cos(#d*)| significantly improves the 
statistical significance for the B° — » D*°7t° mode only. The D*° candidates coming from the decay 
B° D*° 7T° are then required to satisfy | cos(#£>*)l > 0.4 for an efficiency on signal of 90 % and on 
background of 62 % 

An important BB background contribution in the reconstruction of B° — ► D^°tt° comes from 
the color-allowed decay B~ — ► D^°p~. If the charged pion from the decay p~ — > 7r~7r° is not 
used in the construction of the B candidate, D^*'°p~ events mimic the D^°tt° signal. Moreover, 
the BT of B- -> p is almost 50 times larger than the signal BJ 7 , and is known with a 
relative total uncertainty of 13 % for D°p~ and 17 % for D*°p~ [23]. A veto is applied against 
B~ — ► D(*>°p~ events. For each L>(*)°7r candidate, a -B~ candidate is reconstructed in D^ p~ . 
If the B~ candidate mass is in the signal region rriEs(B~) > 5.27 GeV/c 2 , \AE(B~)\ < 100 MeV 
and \m(p~) — ~m(p )pY)Q\ < 250 MeV/c 2 , and if the ir° and Z)° candidates are the same used to 
build the B° candidate, then the B° is rejected. For the reconstruction of -D°7r°, the veto retains 
about 90 % of signal and rejects about 67 % of D°p~ and 44 % of D*°p~ . For the reconstruction 
of D*°7T , the veto efficiency on signal is 80 % while rejecting about 56 and 66 % on D°p~ and 
D*°p~ respectively. 

3.3.4 Choice of one B candidate 

The average number of B candidate per event after all selections ranges between 1 and 1.6 depending 
on the subdecays. The highest multiplicities correspond to modes with the largest number of neutral 
particles in the final state. We keep one B candidate per mode per event. The chosen B is the one 
with the smaller value of: 



'mpVm^Vean' 



°m(D°) / 

+ ( m(h°) - m(fe°)mean \ 2 ^ 



Cm(/i°) 

for D°h modes and 
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' m(D°)-m(D°) mean \ 2 

v a m(D°) J 

' m{h°) - m(/i°)mean \ 2 

v a m(h°) J 

2 



+ / Am- Ammean V ^ 

for the D*°h° modes. The quantities o"m D o and ff mji0 (with respect to m(Z)°) me an and m(/i°) m ean) 
are the resolutions (with respect to means) of the reconstructed mass distributions. The quantities 
Am me an and o& m are respectively the mean and resolution of the Am distribution. These quan- 
tities are obtained from fits of the mass distribution of true simulated candidates selected from 
signal MC simulations. An associated systematic error is calculated later because of the use of MC 
to calculate Ammean and a Am- The efficiency on signal of the choice of one B per event ranges 
from 76 to 100 %. The lower efficiencies correspond to the modes with high neutral multiplicity. 

3.3.5 MC efficiency corrections 

The branching fraction is computed as: 

BT = S -—— 1 (9) 

n BB ' e ' DJ ~ sec 

where BJ- sec is the product of the BJ-'s of the intermediate resonances reconstructed decays, which 
are taken from the PDG [2jj, n B § is the number of BB pairs in data and S is the number of 
signal events remaining after all selections. The quantity e is the signal efficiency of reconstruction 
and selections, computed from the exclusive MC simulations, and the quantity e • BJ-sec is the 
total acceptance. When computing the BT of the sum of Z)W° subdecays, the corresponding total 
acceptances are summed. The total acceptance is corrected by about 1 % (10 %) from the crossfeed 
between D° (D* ) decay modes. The high correction for the sum of D*° decays is mainly due to 
the large crossfeed of D*° — > D°ir° into D*° — > D°j. That correction was obtained by computing 
with exclusive MC's the proportion of intern crossfeed by respect to real signal. 

The selection efficiency from simulation is slightly different from the efficiency in data. The 
MC efficiency for the reconstruction of 7r°/7 is adjusted from the study of r decays in the channels 
r — > p(inr )v T and r — > -kv t |27j . The correction on track efficiency is computed from studies 
of track mis-reconstruction in the decays r — ► (7r + 7r _ )7r _ ^ r and r — > /9 (7r + 7r~)7r~z/ T [28]. The 
simulated efficiency of charged particle identification is compared to the efficiency computed in 
data with pure samples of kaons from the decays D* + — > D°ir + and D° — > K~ir + . The efficiency 
on Kg candidates is modified using a high statistics and a high purity data sample of K® mainly 
arising from the process e + e~ — ► qq. 

The efficiency corrections for the selections applied to Z)W° candidates and on the Fisher 
discriminant for the qq rejection are obtained from the study of the control sample B~ — > D(*^°ir~ . 
That control sample was chosen for its high statistics and purity, and for its similarity with B° — > 
D^°h . The correction is computed from the double ratio eff(data)/eff(MC), where the efficiencies 
are computed from mES fits of B~ — * D(*' 0/ k~ events in data and MC, before and after the 
corresponding selections. The obtained results were checked with the color-allowed control sample 
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B —> Z)W°/9 which has a slightly different kinematic due to the mass of the p and validates 
this correction for modes as D^ rj'. 

3.4 Probability density functions (pdf) and fit procedure 

The signal yield S is extracted by an extended unbinned maximum likelihood (ML) fit of the AE 
distribution in the range -0.280 < AE < 0.280 GeV for m ES > 5.27 GeV/c 2 . Fitting AE allows 
us to model and to fit the complex crossfeed structure without relying on simulation completely. 

Due to the energy loss by photons in the detector material before the EMC, the AE distribution 
for signal is modelled by a modified Novosibirsk function |29j . A Gaussian is added to the modes 
with a large AE resolution to describe the mis-reconstructed events. The signal shape parameters 
are estimated from a ML fit of simulated signal events in exclusive decay modes. 

The dominant crossfeed contribution for B° — ► D°h° comes from the channel B° — ► D*°h° 
when the 7r°/7 from the D*° decay is not reconstructed. Such crossfeed events are shifted in AE 
by about —135 MeV with a tail from D*°(^ D ^f)h in the signal region. 

Table 1: Main crossfeed between signal modes. 



B mode 


Crossfeed mode 




D*°h° 


D* (D ir°)h° 


L>* (L>° 7 )/i , D°h° 


D*°(L> 7 )/i 


D*°(D 7r°)h , D°h° 


D*°h° 


D°h° 



The main crossfeed contributions from the other reconstructed B° modes are shown in Table [TJ 
For each signal mode, different crossfeeds are summed and their contribution is estimated with 
a histogram-based pdf built from the signal MC. The major crossfeed in the reconstruction of 
7T° comes from the decays B — * D^°p (see Section l3.3.3p . Their contribution is 
modelled by a separate histogram-based pdf built from the exclusive MC simulations. 

The combinatorial background from BB and qq are summed and modelled by a 2 nd order 
polynomial. In the cases of B° — ► D^ Q uj(rf) modes, an additional crossfeed contribution comes 
respectively from the decays B~ — > D^ ir^uj(ri). In the reconstruction of B° — ► D°7r° where 
all D° subdecays are summed, an additional crossfeed contribution comes mainly from the decays 
B~ -» D°p~, B° -» Z?°7r vr , B° -» i^Vvr and B° -> D°ir° where D° decays to a CP eigenstate 
or through a semileptonic channel. These peaking backgrounds are significant when summing the 
D° submodes. In these cases, the combinatorial background is modelled by a first-order polynomial 
plus a Gaussian to take into account additional crossfeeds. For B° — > D^ uj(rj) modes, the fraction 
of that Gaussian is left floating in the fit, since the Brs of B~ -> £>W° 7r uj{rj) are not precisely 
known. 

The shape parameters of the combinatorial background pdf are obtained from ML fits on the 
BB and continuum MC, where all signal and crossfeed events have been removed. 

The pdf normalization for the signal, for the combinatorial background and for the crossfeed 
are allowed to float in the fit. The mean of the signal pdf is left floating for the sum of 
subdecays. For each D° submodes the signal mean is fixed to the value obtained with the fit 
on the sum of D° submodes. 
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A given mode B° — » D*°h° can be signal and crossfeed to other modes at the same time. 
In order to use the BJ- computed in this analysis, the yield extraction is performed through an 
iterative fit on successively D*°h° and D°h°. The normalization of crossfeed contribution from 
D(*)°h° is then fixed to the BJ- measured in the previous fit iteration. That iterative method 
converges quickly to a stable value of BJ-, with variation below 10 % of statistical uncertainty, in 
less than 5 iterations. 

We check the absence of bias in our fit by studying embedded toy MC's. The extraction 
procedure is applied to toy samples where background events are generated from the fitted pdf's. 
The toy signal events are taken from the corresponding exclusive MC with a yield corresponding 
to the MC-generated value of the branching fraction BT gen . No significant bias are found. 

3.5 Signal yield extraction 

The iterative fit procedure is applied to data. The fitted AE distributions for the sum of 
submodes are given in Figures [2] and El 




A E (GeV) A E (GeV) A E (GeV) 



Figure 2: Fit of AE distributions in data for modes B° — > D°ir° (top left), B° — > D°uj (top 
middle), B° -> £>°r/(77) ( to P ri ght)> B° L»°r/(7rvr7r ) (bottom left), B° -> D°r]'(inrr]) (bottom 
middle) and B° — > D°i]'(p°'y) (bottom right). The dots with error bars are data, the solid curve 
is the fitted total pdf, the dotted curve is the signal pdf, the dotted-dashed curve is the crossfeed 
pdf, the double dotted-dashed curve is the B~ — > D^*'°p~ pdf and the long dashed curve is the 
combinatorial background pdf. 

The signal and background yields obtained from the fit to data are given in Table [2] with the 
corresponding statistical significances. 
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Figure 3: Fit of AE distributions in data for modes B° -> D*°ir° (top left), B° -> D*°lu (top 
middle), B° -» D* r](Tnnr°) (top right), B° -» D*°rj(~f~f) (bottom left), B° -» D*°(D TT )r] / (p°-f) 
(bottom middle) and B° — ► D*°r]'(inrr]) (bottom right). Detailed legend is provided in Figured 



4 SYSTEMATIC STUDIES 

There are several sources of systematic uncertainty in this analysis. Table [3] summarizes the studied 
systematic uncertainties with their values in percentage of BJ- for the different modes D^*'°h°. The 
large uncertainties for D*°h° mainly come from D* (D°^/)h . 

The categories "vr°/7 detection" ("Tracking") account for the systematics on the reconstruction 
of 7r°/7 and charged particle tracks, and are taken as the uncertainty on the efficiency correction 
computed in the study of r decays (see Section [3. 3. 5p . 

Similarly, the systematic uncertainties on kaon identification and K® selections, are estimated 
from the uncertainties on MC efficiency corrections computed in the study of pure samples of kaons 
and Kg mesons in data respectively (see Section I3.3.5f) . 

The systematic errors on the submodes BJ- is a combination of the uncertainty on each D° 
and h submode [24]. The correlation between the calibration mode D° — > K~ir + and D° — ► X, 
X ^ K~tt + was accounted for. 

The uncertainty related to the number of BB pairs and to the limited MC statistics is also 
included [9]. 

The systematics on resonance mass cut are computed as the relative difference of signal yield 
when the mass mean and resolution for the selection are taken from a fit in data. 

The uncertainties for the qq rejection and the selections are obtained from the study 

performed on the control sample B~ — ► L>(*)°7r~ and are estimated as the uncertainty on the 
efficiency correction eff(data)/eff(MC), including the correlations between the samples before and 
after selections (see Section I3.3.5P . 

The uncertainties for the cuts on p° and D°u> helicities are obtained by varying the selection 
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Table 2: We give the number of signal events (S), crossfeed (-ZV c f), B~ — ► Z)(*)°p~ (iVrj ) and 
combinatorial background (^combi) f rom t ne A-E 1 n t of data, and the i3.F measured for each 
submode. The number of background events iV c f, iVp , and -^comki are computed from the fitted 
pdf in the window \AE\ < 100 MeV. The statistical significance is defined as ^j2ln{Ls+B/ Lb), 
where L$+b is the likelihood of the fit in data for the hypothesis of a signal, while Lb is the 
likelihood for the hypothesis of pure background only. 



B mode 


S 






^combi 


S^xlO" 4 ) 


Statistical 


(decay channel) 










±stat. ± syst. 


significance 


D U 7T U 


3369 ± 102 


135 


440 ± 11 


2585 ± 63 


2.78 ± 0.08 ± 0.20 


35.5 


£>V77) 


1054 ± 49 


44 




823 ± 21 


2.34 ± 0.11 ± 0.17 


26.1 


D v(nmv Q ) 


454 ± 29 


20 




487 ± 15 


2.51 ± 0.16 ± 0.17 


20.3 


D°LU 


1400 ± 58 


68 




1648 ± 28 


2.77 ± 0.13 ± 0.22 


29.4 


D w' (im 77(77)) 


134 ± 14 


6 




74 ± 6 


1.29 ± 0.14 ± 0.09 


14.7 


dV(A) 


290 ± 42 


13 




2382 ± 33 


1.95 ± 0.29 ± 0.30 


7.2 


D*°ir° 


958 ± 73 


121 


1218 ± 41 


844 ± 56 


1.78 ± 0.13 ± 0.23 


15.1 




629 ± 39 


33 




525 ± 18 


2.37 ± 0.15 ± 0.24 


19.4 


D*°r,{TVTT7Y°) 


241 ± 25 


11 




341 ± 13 


2.27 ± 0.23 ± 0.18 


12.2 




1692 ± 86 


60 




4507 ± 47 


4.44 ± 0.23 ± 0.61 


22.3 


D* r)'(irirr)) 


61 ± 10 


3 




34 ± 4 


1.12 ± 0.26 ± 0.27 


8.0 


D*°(D -K )r,'(p°j) 


86 ± 28 


6 




874 ± 20 


1.64 ± 0.53 ± 0.20 


3.3 



cut values by ±10 % around the maximum of statistical significance. 

The category "Signal shape" represents the uncertainty on the shape of the signal pdf for the 
AE fit. The expected shape difference between data and MC is estimated from a study of the 
control sample B~ — > D^*'°p~, hich gives: |Amean| ~ 5.7 MeV and |Acr| ~ 2.3 MeV. The 
systematic on signal shape is obtained by varying the pdf parameters by ±5.7 MeV for the mean 
and ±2.3 MeV for the resolution. 

The uncertainty on the continuum background shape is estimated from the difference of the 
pdf fitted on generic MC's with the pdf fitted in the toes sideband 5.24 < toes < 5.26 GeV/c 2 in 
data. When a Gaussian is added to the combinatorial background shape, the related uncertainty 
is computed by varying its means and resolution by ±lcr. 

The uncertainties related to the branching fractions of the crossfeed from D^*'°h° and D^°p~ 
were computed by varying the BT separately by ±lcr in the fit procedure. The value of a used is 
taken from the fit for D^°h° and from PDG for D^°p~. 

The uncertainty related to the crossfeed shape, modelled by a non-parametric pdf, is obtained 
by shifting and smearing the pdf mean and resolution by ±5.7 and 2.3 MeV respectively, with a 
convolution by a Gaussian. 

The uncertainty in 11 D*°lo helicity" accounts for the hypothesis on the D*°lu polarization in the 
MC generation. As that polarization had never been measured, following HQET and factorization 
based arguments, the D*°lo polarization was simulated using the helicity parameters measured in 
B~ — > D*°p~, thus a longitudinal fraction of Jl = 0.869 [30], confirmed by a CLEO measure- 
ment [31] . That systematic is estimated by a comparison of the efficiency as obtained from an 
exclusive simulation of B° — > D*°lo with a low fi, value of 0.08 and is taken as half of the relative 
difference of the signal yield. 

The most significant systematic uncertainties come from the 7r°/7 reconstruction, the submodes 
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BT, the longitudinal fraction for D*°uj, from the signal and background parametrization on the 
AE fit and from the crossfeed by D^*'°p~. 

Table 3: Systematic uncertainties (%) for the modes B° — * D^°h° where Z)(*) submodes are 
summed. 
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5 RESULTS 

5.1 BT measurement 

The branching fractions obtained from the AE fit on data are given in Table [2] for the sum of the 
£){*)° submodes. The consistency with D° submodes was checked. All the quoted results are the 
branching fraction BTsum computed with the fit of the D° sum. 

The statistical uncertainty c s t a t on BTsum is provided by the fit while the systematic uncer- 
tainty is computed as the average of the systematic weighted by the total corrected acceptance 

_ Si A • Q" S yst,i nn , 
^syst Y~Ai ' ^ > 

where (BTi ± c S y S t,i ± <7 s t a t,j) is the measured branching fraction with the submode i = K~tt + , 
K~7r + 7r~7r + , K~-k + it or K^7r + ir~ . The total uncertainty is eventually computed as the quadrature 
sum of cr s t a t and c S y S t- We checked the compatibility of BT S um with the measurements in the 
submodes BTi for significant signal. The BTs measured with D*° submodes are compatible with 
BT S um, except for B° — > D*°ir° where the incompatibility between D*° — ► D°7 and D*° — > D°ir° is 
due to the large background contribution from D^°p~ . Likewise the measurements in the different 
i] or rl submodes are compatible. The measurements in the rf ) submodes are combined through a 
BT average weighted by the total uncertainty on BT: 
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where (BJ-i±(Ji) is the measured branching fraction with the submode i, and <7j the total uncertainty 
on BT{. The statistical uncertainty on the BT combination is : 



°stat,avg — ^2i/2 ' ( 12 ) 

i I sta.t.i 



1 

'stat,i 

and the systematic is obtained from the Equation (jlOp . 

The ,8 J- J s measured in this analysis are compatible with the previous measurements by BABAR [9\ , 
CLEO [7] and Belle PSHH]. For the modes B° -> L>« ??, and D* (D ir°)Tr°, 

our measured /3.Fs are significantly higher than the ones measured by Belle [TT]. We also 
measured the BT in the dataset studied previously by BABAR [9j and found compatible BJ- values 
with systematic uncertainties lowered by about 25 to 50 %. The improvement mainly comes from 
the selections on £)(*) and on the qq background rejection. 

5.2 Isospin analysis 

The isospin I symmetry relates the amplitudes of the decays B~ — > .D(*) 7r - , 2?° — > D^ + ir~ and 
#o _^ Z) (*)0 7r ) which 

can be written as linear combinations of the isospin eigenstates Ai, I = 1/2, 

3/2 HE3]: 

A(DW+n-) = l/VSA 3/2jDM + ^/3A/ 2 ,DW. 

^pW) = v^sadw-a/iAAiadw, (13) 

and so: 

.A(DW 7r-) = A(DM + ir-) + V2A{D ( * )0 tt°). (14) 

The relative strong phase between the eigenstates A1/2 and A3/2 is noted as 5 for the Dtt system 
and 5* for D*ir system. Final State Interactions between the states D(*)°7t and D(*^ + ir~ may lead 
to a value of 5^ different from zero and, through constructive interference, to a higher value of 
BT for D(*)°vr . One can also define the amplitude ratio R^: 

R< " ~ (15) 

In the HQET limit [361 138] . the factorization model predicts 6^*' = and R^ = 1, while 
SCET [16j [TTJ, [18] predicts 5 to have the same value in the Dtt and D*ir systems. The strong 
phase 5^ can be computed with an isospin analysis of the D^*>ir system. The input parameters 
are the world average values provided by the PDG 2007 partial update [21] of BT{B~ — > D^°tt^), 
BT(W -> DW+7T - ) and of the B lifetime ratio t(B + )/t{B°), the value of BT(B° -> L>(*)°vr°) 
is the one measured in this analysis. We calculate the value of 5^*' and R(*> using a frequentist 
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approach [37] 5 = (30.2j^) and R = 0.70±^ 5 5 for Dtt final states, and 5* = (13.6i? 3 9 6 )°, R* = 
0.641°;°^ for D*tt final states. 

In both D*tt and Dtt cases, the amplitude ratio is significantly different from the factorization 
prediction R(*> = 1. The strong phase is significantly different from zero in the system Dtt, which 
points out that non-factorizable FSI are not negligible. In the system D*tt that phase 5* is smaller 
and compatible with zero due to the lower value of the BJ 7 for Z?*°7r°. 

5.3 Comparison to theoretical predictions on i3jF's 

Table H] compares the BJ- measured with this analysis to the predictions by factorization |131 [21 
[MIES] and pQCD [HI [15]. We confirm the conclusion by the previous BABAR analysis, the values 
measured being higher by a factor of about 4 than the values predicted by factorization. The pQCD 
predictions are closer to experimental values but globally higher except for D°tt°. 

Table 4: Comparison of the measured branching fraction BJ- in units of 10~ 4 , with the predictions 
by factorization and pQCD \14\ I15j. The first quoted uncertainty is statistical and the second is 
systematic. 



B° mode 


This measurement (xlO 4 ) 


factorization (xlO 4 ) 


pQCD (xlO" 4 ) 


L> 7T 


2.78 ± 0.08 ± 0.20 


0.58 ES] 0.70 [2] 


2.3-2.6 




1.78 ± 0.13± 0.23 


0.65 [13] 1.00 [2] 


2.7-2.9 


D° V 


2.41 ± 0.09 ± 0.17 


0.34 Q3] 0.50 [2] 


2.4-3.2 


D*° n 


2.32 ± 0.13 ± 0.22 


0.60 [2] 


2.8-3.8 


D°LU 


2.77 ± 0.13 ± 0.22 


0.66 [13] 0.70 [2] 


5.0-5.6 


D*°uj 


4.44 ± 0.23 ± 0.61 


1.70 H] 


4.9-5.8 


D° V ' 


1.38 ± 0.12 ± 0.22 


0.30-0.32 [35]; 1.70-3.30 [31 


1.7-2.6 


D*°r)' 


1.29 ± 0.23 ± 0.23 


0.41-0.47 [31] 


2.0-3.2 



Factorization predicts the ratio BJ r (B° — > D^°rj')/BJ-(B — > D^°rj) to have a value between 
0.64 and 0.68 [M] linked to the rj — rj' mixing. Table [5] compares that prediction with the experi- 
mental measurements. The measured ratios are smaller than the prediction and are compatible at 
1 a. 

The effective theory SCET [161 El [E] does not predict the absolute value of the BT but that 
the ratios BT(B° -> D*°h°) /BF(B° -> D°h°) ~ 1 for h° = vr°, 77 and rj' . For h° = u that prediction 
holds only for the longitudinal component of B° — > D*°ui as long-distance QCD interactions may 
increase the transverse amplitude. A measurement of fi, is therefor needed. 

The SCET gives also a prediction about the ratio BT{lP -> D^r()IBT(W -> D^°rf) ~ 0.67, 
which is similar to the prediction by factorization. 

The BJ- ratios are given in Table [5j the common systematics between D°h° and D*°h° are 
cancelled. The ratios BT{W -» D*°/i )/jB^(fi° -» D°h°) for h° = vr ,^,^' are compatible with 1. 

6 CONCLUSIONS 

We measured the branching fractions of the color-suppressed decays B° — ► D^°h°, where h° = 
7T , 7], u, and r( with 454x 10 6 BB pairs. Our measurements are in agreement with the previous 
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Table 5: Ratios of branching fractions B^B -> D*°h°)/BJ : '(B The first uncertainty 

is statistical, the second is systematic. 



BT ratio 



This measurement 



BT(D' u tt u ) 

B^(P*°7)( 77 )) 
SJ r (D u »7(77)) 

B^ r (D*°T)(7T7r7r )) 
SJ r (D°77(-!T7r7r u )) 

B^(-D*°r)) 

S.F(.D u u) 
B^t-D* ^' (TTTri))) 

S^(D u r,'(p u 7 ))) 
BJ r (O* T| / ) 

B^(D°r,) 

B.T(D*V) 
8:F(.D*0jj) 



0.64 ± 0.05 ± 0.08 
1.01 ± 0.08 ± 0.08 
0.91 ± 0.11 ± 0.04 
0.96 ± 0.06 ± 0.07 
1.61 ± 0.11 ± 0.20 
0.87 ± 0.22 ± 0.04 
0.84 ± 0.30 ± 0.15 
0.93 ± 0.19 ± 0.11 
0.57 ± 0.06 ± 0.06 
0.55 ± 0.11 ± 0.04 



results [3 O HU] with a significant decrease of both statistical and systematic uncertainties. We 
confirm the significant difference from theoretical predictions by factorization quoted in the previous 
analysis [9] and provide strong constraints on the models of color-suppressed decays. 
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